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between depth within the earth and the velocity of seismic energy. In one embodiment, an model of the relationship between depth 
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A method of processing seismic data 

The present invention relates to a method of processing seismic data such as surface 
seismic data, that is, seismic data acquired using one or more s&xsots disposed on the 
earth^s sur&ce. The sensor(s) may be disposed on land, on the sea-bed, or in the 
land/sea transition zone. The invention may also be applied to processing seismic data 
acquired using a towed sensor array. 

One type of sUr&ce seismic data is seismic data acquired using an Ocean Bottom Cable 
(OBC) seismic surveying method In an OBC seismic survey, a plurality of seismic 
sensors are attached to a support cable at intervals along the length of the cable. The 
cable is then lowered onto the sea-bed to deploy the sensors at their desired locations on 
the sea-bed. OBC seismic sensors generally record the pressure and the elastic 
wavefield of the seismic data. The pressure is a scalar quantity, whereas the elastic 
wavefield is a vector quantity and it is therefore necessary to measure the components 
of the elastic wavefield in three non-coplanar directions, usually in the x-direction, the 
y-direction, and the z-direction (vertical). In total, therefore, four components of the 
seismic data are measured. Multi-component seismic data recording at the sea-floor has 
proven to be a very successful method for imaging through gas saturated overburdens 
and for characterising hydrocarbon reservoirs. 



Another seismic surve)dng arrangement is a towed marine seismic surveying 
arrangement. Such a surveying arrangement uses a seismic sensor, or an army of 
seismic sensors, tiiat is disposed a few tens of metres below the sea surface and that is 
moved through the water by a towing vessel. The seismic sensor, or the array, is 
normally suspended firom a float connected to the towing vessel. A towed marine 
seismic surveying arrangement further includes one or more seismic sources disposed 
below the sea surface, normally suspended jfirom the towing vessel. 

One problem encountered in die processing of sur&ce seismic data.such as OBC 
seismic data is tibat of mode-converted reflections. A mode-converted reflection is a 
reflection in which, for example, a down-going pressure wave (P-wave) is converted on 
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reflection to an up-going shear wave (S-wave). Data obtained from a reflection in 
which a down-going pressure wave (P-wave) is converted on reflection to an up-going 
shear wave are known as P-s data, whereas data obtained from a reflection in which a 
down-going pressure wave is not mode-converted on reflection are known as P-p data. 
The existence of mode-converted reflections (also knows as mode converted events) 
makes it necessary to determine a separate depth-velocity relationship, or velocity field, 
for S-waves, in addition to det^mining a velocity field for P-waves. Other problems 
that arise as a result of mode-conversion reflections include asymmetric paths of seismic 
^ergy, stronger anisotropy effects, and increased attmuation of seismic energy. 

The problem of determining anisotropic velocity fields is of particular importance. The 
anisotropic velocity field used to process the seismic data will affect the moveout 
correction of the data (that is, correcting the data for the horizontal distance - known as 
"offset" - between the source of seismic raergy and the receiver). The anisotropic 
velocity field used to process the seismic data will also affect the shift from midpoiat of 
the coimnon conversion point (CCP). 

One prior art ^proach is to model the velocity field of P-waves and the velocity field of 
S-waves using effective velocity models. In an effective velocity model, variations in 
the propagation velocity of seismic energy that occur within the eartii's interior are 
modelled by an effective velocity of propagation of seismic energy. One well-known 
effective velocity model assumes that tiie medium througji which the seismic energy 
travels is isotropic and that the velocity of seisnuc energy propagating wittiin the 
medium is constant; this is known as Ae *liypeiboIic assumption". Other effective 
velocity models contain non-hyperbolic teems, in an attempt to take account of the 
anisotropy of the earth. 

Even after determining anisotropic effective velocity models for P-waves and S-waves, 
and using these models to moveout correct and stack the P-p and P-s seismic data, there 
remains the problem of coirelating geological events in one set of data with geological 
events in tiieothor set of data. This correlation is required to obtain tiie ratio VpA^s, flie 
ratio of flie votical velocity of P-waves to the vertical velocity of S-waves, which is 
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required for the CCP binning. In turn, the results of (be CCP binning afiect the P-s 
velocity analysis^ so tiiat prior art processing methods require an iterative technique in 
which the results of the CCP binning obtained using an initial velocity analysis are used 
to refine the velocity models re-calculate the ratio VpA^s, and repeat the CCP binning. 
This procedure is time-consuming, and is also highly dependent on the correct 
interpretation of the P-p and P-s data sets in the correlation process. 

A further problrai with the prior art method is that in some cases reliable velocity fields 
cannot be determined. This can occur in gas clouds, for example, or if there is severe 
contamination of the seismic data by multiple reflections. Moreover, even if reliable 
effective, anisotropic velocity fields can be detennined for both P-waves and S-waves 
fi-om towed marine seismic data or surface seismic data such as OBC seismic data, it is 
not possible to determine a single, consistent anisotropic effective velocity model firom 
surface data or towed marine seismic data alone. 

A further problem associated with processing surface seismic data or towed marine 
seismic data is that mode conversion is not limited to occurring upon reflection. Mode 
conversion may also occur when, for example, a down-going wave passes firom one 
layer of the earth's interior to anothCT layer having different properties, and in particxilar 
when a down-going wave passes into a layer in which the velocity of seismic energy is 
high. Thus a down-going P-wave may be mode-converted to a down-going S-wave at a 
layer boundary within the earfli rathar than on reflection. A down-going P-wave that is 
mode-converted to a down-going S-wave and is subsequratly reflected to give an iqp- 
going S(?)-wave is referred to as a P-S-s mode or P-S-s event The presence of P-S-s 
modes complicates the interpretation of seismic data relating to a target reflector that is 
disposed below a high velocity layer. 

One prior art approadi to imaging a target reflector beneath a high velocity layer has 
considered coupling between P-waves and S-waves in the high velocity layer and pure 
P-wave propagation below the high velocity layer - see, for example, G.W.Pumell, 
"Geophysics'', Vol 57, No 1 1 (1992), ppl444.1452. An alternative approach has 
considered a P*s reflection in which a down-going P-wave is mode-converted to an up- 
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going S-wave on reflection - see, for example RJLKendall ct al in Expanded Abstracts 
of 68* Annual Ihtemational Meeting of Society of Exploration Geophysicists (1998), 
pp2052-2055- However, there is currently no satisfactory approach for dealing with 
P-S-s mode converted events. 

A first aspect of the present inv^tion provides a method of processing surface seismic 
data or towed marine seismic data, flie method comprising the step of using vertical 
seismic profile (VSP) seismic data to determine a model of flie velocity of seismic 
energy. A second aspect of the invention provides a method of processing surface 
seismic data, the method conqirising the step of using vertical seismic profile (VSP) 
seismic data to determine a model of tihe velocity of seismic energy. The invention 
enables a more accurate model of tiie relationship between the depth and the velocity of 
seismic energy to be detennined. 

In a preferred raibodiment, the method comprises the step of calibrating a velocity 
model for the surface seismic data using VSP seismic data. 

In this embodiment an initial velocity model is calibrated using VSP seismic data. The 
initial velocity model may be derived, for example, firom seismic data for a dipole sonic 
log. Calibrating the initial velocity model using VSP seismic data enables an accinate 
velocity model to be produced at the start of the analysis process. This shortens, or 
even eliminates, the iterative process described above. Moreover, the invention also 
allows a reliable velocity model to be detennined even if the suifece seismic data are 
contaminated by multiple reflections or by gas clouds. 

In a preferred embodiment, the method comprises the further step of calibrating the 
velocity model using surfece seismic data. This finther calibration step may be used, 
for exan^le, to provide increased accuracy of the velocity model at dep^ above the 
shallowest VSP sensor. Alternatively, towed marine seismic data may be used in this 
calibration step. 
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In a prefexied embodiment, a first velocity model for P-waves and a second velocity 
model for S-waves are determined using the VSP seismic data. Thus, the present 
invention enables accurate processing of seismic data which contains mode conversion 
events. 



A further embodiment of the invention comprises the step of processing the seismic data 
using the or each velocity model detennined using the VSP seismic data. 

The invention may be applied to OBC seismic data. However, it is not limited to 
plication to OBC seismic data, but may be applied to all surface seismic data or to 
towed marine seismic data. It may be ^plied to multi-component seismic data or to 
single component seismic data 

Other preferred features of the present mvention are set out in the dependent claims. 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying Figures in which: 

Figure 1 is a schematic illustratioa of a seismic survey; 

Figure 2 is a schematic Qhistration of a velocity-dq)th model obtained using a method 
of die present invention; 

Figure 3 shows data traces simulated using the model of Figure 2; 

Figure 4 shows a comparison between the results obtained by processing seismic data 
using an isotropic velocity-depth model and the results obtained using a velocity-depth 
model of the present invention; and 

Figure 5 shows the results of processing seismic data using a velocity-dq)th model 
according to the present invCTtion. 
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Figure 1 is a schematic view of an OBC seismic surveying axrangement The 
arrangement has a seismic source 1 and a seismic sensor 2 (otherwise known as a 
receiver) which togetfier form a surface seismic survey. In the embodunent shown in 
Figure 1 the seismic sensor 2 is a multi-compon^t OBC sensor that forms part of an 
OBC sensor array» and is disposed on the sea*bed 3. Other sensors in the OBC sensor 
anray, and the support cable for the OBC sensors, are omitted jfrom Figure 1 for clarity. 
As noted above, however, the invention is not, however, limited to processing data 
acquired in an OBC seismic survey. 

The seismic source 1 and the OBC sensor 2 are being used to obtain data about a target 
reflector 4 disposed below the sea-bed. The path of the primary reflection of seismic 
©tiergy from the seismic source 1 to the OBC sensor 2 via the reflector 4 is indicated in 
Figure 1 as the path 5. 

The geological structure of the earth above the target reflector 4 of Figure 1 is not 
uniform. The earth's interior between the sea-bed 3 and the target reflector 4 includes 
an layer 6 in which the velocity of seismic energy is significantly greater than in the 
layer 7 below the high velocity layer 6 or in the layer 8 below the high velocity layer 6 
(only three layers 6,7,8 are shown between the sea-bed and the target reflector 4 in 
Figure 1, but in practice fliere will be many differ^t geological layers between the sea- 
bed and die target reflector 4). The high velocity layer 6 might be, for example, a chalk 
layer, possibly a weathered chalk layer. One effect of the high velocity layer 6 is that 
down-going seismic energy undergoes mode conversion at the upper surface of the 
high-velocity layer 6; A down-going P-wave undergoes mode conversion to a down- 
going S-wave at the upper surface of the high-velocity layer 6, and is then reflected at 
the target reflector to produce an up-gomg S-wave. Qn practice some reflection and 
refraction will occur at the layer bo\mdaries, and seismic energy paths which involve 
multiple reflections may weU exist Reflections and refraction at the layer boundaries 
have been omitted from Figure L) 



As explained above conventional methods of processing the smiace seismic data 
acquired by flie sensor 2 would require a time-consuming it^tive technique^ even in 
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the absence of mode conversion at the upper surface of the high- velocity layer 6. For 
die geological structure shown in Figure 1, the processing of seismic data relating to the 
target reflector 4 will be complicated by the P-S-s mode conversions that arise due to 
the high velocity layer 6. 

According to the processing method of the present invention, data from a vertical 
seismic profile (VSP) seismic survey is used in processing the data obtained by the 
surfece sensor 2. More specifically, in an embodiment of the invention VSP data is 
used to detemune a model of the relationship between depth and the velocity of seismic 
energy, and this velocity model is then used in the processing of tte seismic data 
acquired by die sensor 2. 

The surveying arrangement for obtaining VSP seismic data is also shown in Figure 1. 
The VSP surveying arrangement consists essentiaUy of a sdsmic source 9 disposed just 
under the surface of the sea (generally suspended fi'om a survey vessel) and seismic 
sensors 10 disposed at different depths within a borehole 1 1. Seismic energy emitted by 
the VSP seismic source 9 may travel to one of the sensors 10 either by a direct path 
througji the overburden (indicated as path 12 in Figure 1) or via a reflection at the target 
reflector 4 (as indicated by path 13 in Figure 1). By "overburden** is meant the part of 
the earth's interior that extends &om the sea-bed 3 down to the depth of the shallowest 
of the VSP seismic sensors 10. The effects of refi:action and reflection at the int^&ce 8 
have been omitted from paths 12 and 13. 

Three seismic sensors 10 are shown in the borehole 1 1 of Figure 1, disposed at three 
dififerent depths, but the invention is not limited to this VSP survey arrangement 

Figure 1 shows an offset VSP seismic survey in which there is a fixed, non-zero oSsdt 
between the source 9 and the sensors 10, but it is also possible to carry out a VSP 
survey with an o&et of ^proximately zero. It is also possible to cany out a VSP 
survey at two or more different ofi&ets; a VSP survey in which data is available for 
more ttian one value of ofEset is generally known as a 'Valkaway" VSP seismic survey. 
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The gmeral tenn •*VSP seismic data" as used herein is intended to include zero-ofifeet 
VSP data, fixed offset VSP data and multi-offset (walkaway) VSP data. 

The initial step in the processing method of the present inv^tion is to develop an initial 
model of the relation between depth and the velocity of seismic energy for P-waves. 
This is done, in this embodiment, by extending data from a dipole sonic log through the 
overburden. This can be done using one or more of a compaction trend velocity, 
regional check-shot information, or surface seismic Dix velocities. Once the extended 
model for the velocity of P-waves has been established it is then calibrated using a 
standard process of sonic cahbration (that is, using direct VSP arrival times obtained 
from zero-of&et or low-o£&et VSP seismic data). 

Once a satisfactory model for the velocity of seismic energy has been developed for P- 
waves, the next step is to detemiine a model for the velocity of seismic energy for S- 
waves. The initial velocity model used for S-waves is obtained from the shear sonic, 
and this is calibrated using S-wave velocities picked from the horizontal components of 
zero-oflfeet VSP data or ofeet VSP data. 

It has been found that in the case of a vertical borehole there is little difiference between 
the S-wave depth-velocity relationship derived from a dipole sonic log and the S-wave 
depth- velocity relationship derived from VSP seismic data. This is not the case, 
however, for a borehole that is deviated from the vertical. For a borehole that is 
deviated from the vertical, polar anisotropy can cause the S-wave velocity derived from 
a dipole sonic log to be significantly higher than the true S-wave velocity. In flus case, 
the value for the S-wave velocity derived from the VSP data is used instead of the S- 
wave velocity doived Scorn the dipole sonic log. 

No reliable direct measurement of the S-wave velocity for depths above the shallowest 
VSP seasoT 10 is available, since it is not always possible to identify amode conversion 
from P-waves to S-waves at the sea-bed from VSP data. The S-wave velocity for 
depths above the shallowest VSP sensor 10 may be estimated Scorn the P-wave velocity 
measured for the oveiburden, using a known or estimated relationship between tfie P- 
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wave velocity and the S-wave velocity. In a preferred embodiment, however, the 
relationship between depth and S-wave velocity for depths above the shallowest VSP 
sensor 10 is detennined using an extrapolated compaction trend in VpA^s. The 
compaction trend in VpA/s is constrained to pass through a predetennined value at the 
sea-bed 3, and varying the value of VpA^s provides a method caUbrating the vertical S- 
wave velocity for depths above the shallowest VSP sensor 10. This calibration uses the 
near offset moveout of shallow P-s mode conversion reflections in OBC radial data; 
ray-based moveout correction is used to remove the near offeet reflection moveout by 
varying the VpA^s ratio and thus varying the S-wave velocity m the oveiburden. This 
quickly leads to a model for the S-wave velocity in the overbunl^ 

The resulting model for the vertical S~wave velocity is preferably then conqyared with 
P-s mode converted events in of&et VSP data, and adjusted if necessary. 

Preferably, the method includes a further calibration step of the S-wave velocity for 
dq>ths above the shallowest VSP sensor 10. In this liirther cahbration step, shallow 
marker events on the corrected OBC common receiver gather are compared with VSP 
data for P-s mode converted reflections, or with simulated data for P-s mode converted 
reflections. The model for the S-wave velocity for the overburdm is further adjusted, if 
necessary, to improve the results of this comparison. In order to carry out this step by 
comparing the OBC data with simulated data for P-s mode conv^ed reflections a 
density model of the geological stracture in the area of the seismic survey is required, 
and this can be developed fiom the edited density log and jfrom a calibrated Gardner's 
relationship (a calibrated eo^irical relationship between the P-wave velocity and ^e 
density) for the overburden. 

Thus, the present invention makes it possible to develop a velocity model for S-waves, 
by using VSP data and surface seismic data (in this embodimmt OBC seismic data) to 
calibrate an initial model developed fiom borehole seismic data. The velocity model 
can then be used to process the surface seismic data 
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In a further embodimoit of the present invention, the velocity model is further refined to 
take account of polar anisotropy. It is possible to detennine polar anisotropy by 
min i mi sing the residual long-offset moveout of P-p and P-s reflections on OBC data. 
Alternatively, if data from a walkaway VSP seismic survey is available then polar 
anisotropy may be determined by using the arrival times of the direct seismic pulses 
recorded in the walkaway VSP seismic survey. 

The step of calibrating the velocity model for polar anisotropy using walkaway VSP 
data is carried out using a simple two parameter time inversion. The ellipticity 
parameter and the anellipticity parameter are each assumed to vary linearly with depth 
from a value of zero at the sea-bed. ITie ellipticity and anellipticity paiamet^ are 
varied, and the values that minimise the diflference between the modelled direct arrival 
time and ^ measured direct anival time of P-waves are selected. 

The above description has treated the overburden as a single layer. The ov^urden 
will, however, g«ierally consist of more than layer, and the layers fliat extend 
throughout the well dq>th are each calibrated using VSP seismic data (sonic 
calibration). Each layer that extends through the well dqpth is also preferably calibrated 
for anisotropy using moveout times, using apiecewise gradient assmnption, to obtain 
the intrinsic VTI parameters for the layer. 

It is preferably assumed that a layer is isotropic if the VpA^s ratio for that layer is below 
a threshold vahie. This enables accurate modelling of a layer that is isotropic, or only 
weakly anisotropic, for example such as a reservoir sand layer. 

The direct travel times are preferably corrected for dip, for example using the method 
proposed by CM. Sayers in "Geophysics" Vol 62. No 3 (1997), pp723-729. 

The layered VTI model developed according to the above-described method of the 
invention is valid for convated Sv modes because, provided that the vertical velocities 
of P-waves and S-waves have been correctly determmed, long ofifeet moveout is 
sufficieaQt to estimate the VTI parameters governing the propagation of Sv modes. If 
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VSP walkaway data is available for receivers at more than one depth, a piecewise 
gradient model in anisotropy parameters can be detemiined 

In a further embodiment of the invention^ the velocity model is fiuther calibrated to take 
account of anelastic attenuation, generally know as Q (a dimensionless parameter). This 
calibration can be done by inverting direct arrival wavefoims obtained using zero-oflfset 
VSP and/or walkaway VSP, for example using the coherency Q inversion method 
described by W.S. Leaney in 'Walkaway Q Inversion" in Expanded Abstracts of 69*^ 
Annual International Meeting of Society of £;q)loration Geophysicists (1999), pp6-52. 

In one study it was found that a single Q value, of 54, provided good compensation of 
VSP waveforms at all depths, although the chalk interval indicated a slightly lower Q 
vahie for the chalk layer than for ottier layers. If offset VSP data are available with 
sufficient ofifset to generate significant mode conversions, then a reliable estimate of Q 
for S'waves, (^hear, is also possible over the depth interval of ttie VSP data. 

Figure 2 illustrates an anisotropic velocity model obtained by a method of the present 
invention. This model has 160 layers, and extends to a depth of 2,900 metres. Starting 
Scorn the left. Figure 1 illustrates flie vertical velocity of S-waves, Vs; the ratio between 
the vertical velocity of P-waves and the horizontal velocity of P-waves Vp(v)A^p(h); the 
ratio of the vertical velocity of P-waves to the vertical velocity of S-waves, Vp/Vs; 
Schoenberg's anisotropy parameter, and the density. The model uses a constant Q 
value of 54 for both P-waves and S-waves. 

The 160 layer model shown in Figure 2 was used to simulate data, although a simpler 
model having only 25 layers was used in processing to compensate for Q and geometric 
spreading, for nomoal moveout ^^MO) correction, for CCP binning and for angle band 
mute design. 

The ^proach described above based on a VTI velocity model produced good results in 
processing P-s mode conversion events. The model may also be used to process P-S-s 
mode conversion events. Simulations using the model of Figure 2 showed diat flie most 
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efficient mode conversion from a down-going P-wave to a down-going S-wave 
occurred at the upper surface of a high velocity layer (such as a chalk layer). This 
conversion mode is advant^eous, since the lower velocity once the mode conversion 
has occurred should lead to better vertical resolution and to less contamination of the 
data by multiple reflections. 

Figure 3 shows data traces simulated using the model of Figure 2. Figure 3 shows, from 
the left, the vertical velocity of S-waves, the vertical velocity of P-waves, and the 
simulated seismic traces (showing only primary reflections). The y-axis of Figure 3 
r^resents two-way travel time (in seconds). The x-axis is in metres/second for the 
velocity of S-waves and P-waves, and in metres for 4e sunulated seismic traces. A 
high velocity chalk layer exists, with the upper surface of the high velocity layer 
corresponding to a travel time of 3.12s. Above the upper surface of the chalk layer the 
events are P-s reflections, and below the upper siirface of the chalk layer the events are 
P-S-s reflections. The data was simulated using a 10-45® angle band for P-s stacking 
and a 20** inner mute for P-S-s stacking to avoid the polarity reversal that occurs on 
reflection of a down-going S-wave. 

Figure 4 illustrates the effect of anisotropy on P-S-s reflection evmts occurring below a 
high velocity layer. The x-axis rqiresrats horizontal distance (in metres), with the 
borehole located at x*0. The y-axis represents P-s two-way times (in seconds). A 
hig^h velocity chalk 1^^ exists, with the upper surface of the chalk layer corresponding 
to a travel time of 3.12 seconds. 

On fbs left of the borehole the Figure shows the results of moveout correcting P-S-s 
events using an isotropic model, and on the rigjit of the borehole the Figure shows the 
results of moveont correcting P-S-s evwts using an borehole-calibrated Vn model of 
the present invntion. The borehole-calibrated VTl model of the invention incorporates 
geometric spreading, and also has Q compensation (using a Q value of Q - 54 for botti 
S-waves and P-waves). Thereadualmoveout,relative to thebordioleK^aUbratedVTI 
model of the present invention is indicated on the Figure, and it will be seen that the 
residual moveout for the isotropic model inoeases strongly with increasing ofiset 
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Once the surface seismic data has been moveout corrected using a velocity model of the 
present invention, common reflection point binnning (CRP binning) is carried out This 
is generally analogous to CCP binning for a P-s event. Other processing operations, 
such as deconvolution and correction for receiver statics may also be apphed. 

Figure S shows the result obtained from a velocity model of the present invention. The 
horizontal scale represents the common conversion point, and the vertical scale 
represents depth in km. P-s mode events caused by reflections at the upper surface of a 
high velocity chalk layer are prominent at around 1.7km depth, but events caused by 
reflection at the base of the chalk layer, or by reflection below the chalk layer, are 
clearly visible. The section below the chalk layer shows better vertical resolution than 
the section above the chalk layer. This Figure clearly shows tiiat the present invention 
allows high-quality imaging below a high velocity layer. 

Figure 4 also shows (as the white trace) the relative shear impedance log at the location 
of the borehole. 

The embodiments described above relate to a caUbrating a velocity model using a one- 
dimensional VTI velocity model. This model is strictly valid only at locations near the 
borehole. There are a niunber of ways in which the model could be extended to 
distances further away from the borehole. In one spproach, ihe VTI velocity-depth 
model could be averaged for effective time-based parameters to seed parameter scans 
away from well locations, using the method proposed by I. Tsvankin et al in 
"Geophysics" Vol 519 (1994) ppl290-1304, or using the method proposed by T. CSieret 
et al in a p^er submitted to the 70^ Annual Intmiational Meeting of Society of 
Exploration Geophysicists. 

Alternatively, if additional well data are available in the volume of tiie earth covered by 
the sur&ce seismic data, the method described above can be x&pestcd at additional 
control points, and a locally one-dim^isional model can be interposed between ttie 
control points. 
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A further approach to extending the model would be to use flie one-dimensional model 
to constrain a 3-dimeasional anisotropic prestack depth migration process. 

The embodiment of the invention described above relates to the use of VSP data in the 
processing of surface seismic data, for example OBC seismic data^ that contains mode 
conversion events. The invention is not limited to this specific application, however, 
but may be ^plied to all sur&ce seismic data or towed marine seismic data. The 
invention may be applied to multi-component seismic data or to single component 
seismic data. 

A further application of the invention is to the long-offset processing of surface seismic 
data. Processing of long-oJBset seismic data requires accurate knowledge of the 
anisotropy parameters, of tiie effective value of the anelasticity parameter Q, and of the 
AVO response. This information may be obtained in a multi-oflfeet VSP seismic survey 
using a three-component seismic receiver located within a borehole, and used to adjust 
an initial velocity mode to improve its accuracy. The data may then be processed using 
the iniproved velocity model. 

A further ^plication of the invention is to pre-stack depth migration of seismic data. 
Accurate migration requires an accurate velocity model, which conventional techniques 
are not always able to provide. Once an initial velocity model has been obtained, the 
method of the invention may be used to improve the accuracy of this velocity model, by 
constraining the velocity model so that it rq>roduces travel times obtained fix>m VSP 
data. This improves the accuracy of the velocity model, particularly in the vicinity of 
the borehole used to obtain the VSP data. Once the velocity niodel has been corrected 
in this way, it can be further corrected by, for example;, including anisotropy in the 
model [does the anisotropy information also come from VSP data?]. This method is 
particularly effective if 3D VSP data is used to constrain the initial velocity model. 
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CXAIMS: 



1 . A method of processing surface seismic data or towed marine seismic data, the 
method comprising the step of using vertical seismic profile (VSP) seismic data to 
determine a model of the velocity of seismic enei^. 

2. A method of processing surface seismic data, the method comprising the step of 
using vertical seismic profile (VSP) seismic data to determine a model of the velocity of 
seismic en^gy. 

3. A method as claimed in claim 1 or 2 and comprising the step of calibrating a 
model of the velocity of seismic energy using VSP seismic data. 

4. A method as claimed in claim 3 and comprising the further step of calibrating 
the model of the velocity of seismic energy using surface seismic data. 

5. A method as claimed in claim 1, 2, 3 or 4, wherein the method comprises 
determining a first model of the velocity of seismic energy for P-waves and determining 
a second model of the velocity of seismic energy for S-waves. 

6. A method as claimed in any preceding claim wherein the or each velocity model 
is an anisotropic velocity model. 

7. A method as claimed in claim 6 wherein the or each velocity model is a VTI 
velocity model. 

8. A method as claimed in any preceding claim wherein the or each velocity model 
is a one-dimensional velocity model. 

9 . A method as claimed in any preceding claim wherein the step of calibrating the 
or each velocity model using VSP seismic data conqxrises calibrating the or each 
velocity model for the effects of anelastic attenuation of seismic energy. 
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10. A method as claimed in any preceding claim and further comprising the step of 
processing the seismic data using the or each model of tibe velocity of seismic energy 
detemuned using the VSP seismic data. 

11. A method as claimed in any preceding claim wherein the seismic data is Ocean 
Bottom Cable (06C) seismic data. 

12. A method as claimed in any preceding claim wherein the seismic data is multi- 
component seismic data. 



